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O 

■ We discuss the SU{2) chiral sigma model in the context of nuclear matter using a mean 
' field approach at high density. In this model we include a dynamically generated isoscalar 

vector field and higher-order terms in the scalar field. With the inclusion of these, we 
reproduce the empirical values of the nuclear matter saturation density, binding energy, and 
nuclear incompressibility. The value of the incompressibility is chosen according to recently 
obtained heavy-ion collision data. We then apply the same dynamical model to neutron- 
rich matter in beta equilibrium, related to neutron star structure. The maximum mass and 
(N , corresponding radius of stable non-rotating neutron stars are found to be in the observational 

limit. 

(N 

> . 

OO , §1. Introduction 

■ It has been argued that chiral symmetry is a good hadron symmetry^) which 
i ranks only below isotopic spin symmetry. The spontaneous breaking of this sym- 
' metry generates the constituent quark masses and hence various hadron masses, 

including the nucleon mass. Therefore, the theories of dense nuclear matter, where 
dynamical modifications of hadron masses are expected should possess this symme- 
try. In recent decades, three-body forces in the equation of state at high density have 
^ ■ been studied by several authors. These studies also support the importance of 

pH I studying the chiral sigma model, because the non-linear terms in the chiral sigma La- 

grangian can give rise to three-body forces. Theories of dense nuclear matter should 
be capable of describing the bulk properties of nuclear matter, such as the binding 
energy per nucleon, saturation density, compression modulus/incompressibility and 
' symmetry energy. Presently there is no theory of dense nuclear matter that describes 

all the nuclear matter properties and possesses chiral symmetry. 

A chiral Lagrangian using a scalar (sigma) field was originally introduced by 
Gell-Mann and Levy, '^^ and later the importance of chiral symmetry in nuclear mat- 
ter was emphasized by Lee and Wick. The usual theory of pions does not possess 
the empirically desirable saturation properties for nuclear matter. For this reason, 
an isoscalar vector field with a dynamically generated mass was introduced into 
the theory of nuclear matter and it enabled us to have a saturation density in nu- 
clear matter equation of state. In the standard sigma model, the value of the 
incompressibility parameter of nuclear matter turns out to be quite large — several 
times the desirable value — and can be reduced only by introducing the scalar field 
self-interactions with adjustable coefficients. Several papers ^^''^^ have attempted to 



X 



JSPS fellow. 



typeset using pTpJ^sty <ver.l.O> 



2 Pradip Kumar Sahu and Akira Ohnishi 

derive the chiral sigma model equation of state at high matter density with normal 
nuclear matter saturation as well as a desirable incompressibility value. However, 
in these theories, the mass of the isoscalar vector field is not generated dynamically. 
This fact can be considered a shortcoming of the chiral symmetry model. 

A few years ago, we proposed the SU{2) chiral sigma model to describe the 
properties of nuclear matter. In that work we adopted an approach in which the 
mass of the isoscalar vector field is generated dynamically. To ensure the saturation 
properties of nuclear matter, inclusion of such a vector field is necessary. The nucleon 
effective mass thus acquires a self-consistent density dependence both on the scalar 
and the vector meson fields, and these meson fields are treated in the mean-field 
theory. To describe the nuclear matter properties, we have two parameters in the 
theory: the ratio of the coupling constants to the scalar and to the isoscalar vector 
fields. In this approach the value of incompressibility at saturation density is rela- 
tively large, which is an undesirable feature, as far as nuclear matter at saturation 
and higher densities is concerned. In the present calculation, we rectify the above 
mentioned shortcoming at the nuclear matter saturation density by including higher- 
order terms of the scalar field potential in our proposed chiral sigma Lagrangian. ^) 
In this way, we get two extra parameters in the mean-field approach. These are 
fitted with the nuclear matter properties at the saturation density. 

The paper is organized as follows: In §2 we briefly describe the proposed SU{2) 
chiral sigma model and the derivation of the equation of state for the nuclear matter 
density. Section 3 contains the neutron star matter equation of state in beta equi- 
librium with the neutron star results. A conclusion and summary are presented in 
§4. 

§2. The SU{2) chiral sigma model and equation of state 

The Lagrangian for the SU{2) chiral sigma model can be written as (we choose 
units in which h = c = 1) 

+ ^gjx^u;^ij^ + j^^gJiio^uj^f - Da, (1) 

where Ff^^ = d^j,u}y — di^io^ and a;^ = tt^ + o"^ , ^ is the nucleon isospin doublet, tt is 
the pseudoscalar-isovector pion field, a is the scalar field, and D is a constant. The 
Lagrangian includes a dynamically generated isoscalar vector field, a;^, that couples 
to the conserved baryonic current = tp^f^ip- B and C arc constant parameters 
included the higher-order self-interaction of the scalar field in the potential. In the 
fourth order term in the omega fields, the quantity ^ is a constant parameter. For 
simplicity, we set ^ to zero in our calculation. 



SU (2) Chiral Sigma Model and Properties of Neutron Stars 



3 



The interactions of the scalar and the pscudoscalar mesons with the vector me- 
son generate a mass of the latter through the spontaneous breaking of the chiral 
symmetry. The masses of the nucleon, the scalar meson and the vector meson are 
respectively given by 

M = gaXo, rua = V2Xxo, = g^^Xo , (2) 

where Xo is the vacuum expectation value of the a field, A = (mg-^ — m^^)/(2/^^), 
with TTT-TT the pion mass and the pion decay coupling constant, and g,^ and g^ are 
the coupling constants for the vector and scalar fields, respectively. Throughout this 
paper, we limit ourselves to a mean field treatment and ignore the explicit role of tt 
mesons. 

The equation of motion of fields is obtained by adopting the mean-field approxi- 
mation. This approach has been used extensively to obtain field theoretical equation 
of state models for high density matter. Using the mean-field ansatz, the equation 
of motion for the isoscalar vector field is 

Wo = 2 ' (^) 

QojX^ 

and the equation of motion for the scalar field in terms oi y = xjxg is 



2c^c^ _ r'^F Pdk ^ 
where M* = yM is the effective mass of the nucleon and 

— 2/2 _ 2 / 2 /r\ 

The quantity ub is the baryon density, or equivalently, ub = lo^ d^^i where 
is the Fermi momentum and 7 is the nucleon spin degeneracy factor. 
The equation of state is calculated from the diagonal components of the con- 
served total stress tensor corresponding to the Lagrangian together with the mean- 
field equation of motion for the fermion field and a mean-field approximation for 
the meson fields. The total energy density, e, and pressure, P, of the many-nucleon 
system are the following: 



27r^ Jo 
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The energy per nucleon is E/A = e/ns, 
where 7 = 4 for symmetric nuclear mat- 
ter. 

In the above equations, we have four 
parameters: the nucleon coupling to the 
scalar and the vector fields, and c^, 
and the coefficients in the scalar poten- 
tial terms, B and C. These are ob- 
tained by fitting the saturation values of 
binding energy/nucleon (—16.3 MeV), 
baryon density (0.153 fm~^), and effec- 
tive (Landau) mass (0.85M). The nu- 
clear incompressibility is somewhat un- 
certain at saturation and therefore we 
choose it between 250 MeV and 350 
MeV, i.e. ^ 300 MeV, in accordance 
with recent heavy-ion collision data. 
11), 12) These values are = 1.9989 fm^, 
= 6.8158 fm^, S = -100 and C = 
— 133.6. In Fig. 1, we present the en- 
ergy/nucleon as a function of the baryon 
density. The solid curve (MCH) corresponds to the modified chiral sigma model pre- 
sented above, where the nuclear incompressibility is around 300 MeV. This equation 
of state is much softer than that of the original chiral sigma model (CH), which is 
represented by the dotted curve. This is due to the additional terms B and C in the 
modified chiral sigma model in the potential term in equation (1). For comparison, 
we display the dashed curve (NL3), which was derived from recent heavy- ion colli- 
sion data with incompressibility around 340 MeV. We notice from Fig. 1 that the 
MCH model is softer than the NL3 model at high density. This is due to the slight 
difference in the incompressibility. 

§3. Equation of state in neutron star matter 

In the interior of neutron stars, i.e. at high density, the neutron chemical poten- 
tial exceeds the combined mass of the proton and electron. Therefore, asymmetric 
matter, with an admixture of electrons, rather than pure neutron matter, is a more 
likely composition of matter in neutron star interiors. The concentrations of neu- 
trons, protons and electrons can be determined from the condition of beta equilib- 
rium (n -ir^ p + e + i/) and from charge neutrality, assuming that neutrinos are not 
degenerate. We have 

IJ'n = IJ-P + fJ'e, rip = He, (7) 

where /x^ is the chemical potential of particle species i. For the purpose of describing 
neutron-rich matter, we include the interaction due to the isospin triplet p meson in 
the Lagrangian (1). The following terms are included in the Lagrangian: 




-50 

1 2 3 4 5 

Fig. 1. Energy per nucleon as a function of 
baryon density in units of no. The solid 
curve (MCH) corresponds to modified chi- 
ral sigma model, where the nuclear incom- 
pressibility is around 300 MoV, the dashed 
curve (NL3) corresponds to the equation of 
state derived from recent heavy-ion colli- 
sion data with incompressibility around 
340 MeV, and the dotted curve corresponds 
to the original chiral sigma model, with 
very high incompressibility. 
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-^G,, ■ G^^ + ^m>^ • - ^gp^ip^ ■ r^)^ . (8) 
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Here G^u = d^^^Py — d^p^. Using the mean-field approximation in the equation of 
motion for p, the fohowing density dependence equation is obtained: 



„3_ 9p 

-I 

Prom the semi-empirical nuclear mass formula, the symmetric energy coefficient is 



Pl = 4^^{n,-nn). (9) 



where Cp = (/p/j^n^ and fci;' = (67r^nB/7)^^^(ras = + We fix the coupling 
constant Cp by requiring that agym correspond to the empirical value, 32 MeV. This 
gives Cp = 4.66 fm^. The inclusion of the p meson in the Lagrangian will contribute 
the term m?p{p'^)'^ /2 to the energy density and pressure. Then the equation of state 
for neutron rich nuclear matter in beta equilibrium is calculated using the expression 
for the energy density e and the pressure P as follows: 



M2(1- 7/2)2 5 



2\3 I Ci^n^ 



In these equations epc a-^d Pfg a-^'e the relativistic non- interacting energy density 
and pressure of the baryons (with effective masses) and electrons (i), respectively. 

Figure 2, displays the pressure versus the total mass-energy density for neutron- 
rich matter in beta equilibrium. The solid curve corresponds to the modified chiral 
sigma model MCH considered in the present calculation. Inclusion of the two param- 
eters, B and C in the potential term, gives a reasonable value for the incompressibility 
at the saturation density of nuclear matter. Hence, the MCH model is much softer 
than the original CH model with regard to the neutron star matter equation of state. 
The physics behind this is that the pressure generated by the self-interaction of scalar 
fields at high density is less; i.e., the pressure with density decreases more than that 
for the CH model. The dashed curve (NL3) represents the neutron rich equation of 
state, which was derived from the heavy-ion collision data. From Fig. 2, we note, 
in comparison with the NL3 model, that the MCH model is stiffer in the low density 
range and becomes softer above three times the normal nuclear matter density. 
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The mass and radius of a neutron star are characterized by its structure. These 
are determined from the equations that describe the hydrostatic equihbrium of degen- 
erate stars without rotation in general relativity, the Tolman-Oppenheimer-Volkoff 
(TOV) equations: ^4)- 



dp G{e + p/c^){m + 47rr^p/c^) dm 
dr r^{l — 2Gm/rc^) ' dr 



= 47rr^e 



(12) 
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o 
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Fig. 2. The neutron star matter pressure as a 
function of the energy density. The models 
are the same as in Fig. 1. 



Here p and e are the pressure and to- 
tal mass-energy density, and m(r) is 
the mass contained in a volume of ra- 
dius r. The quantity G is the gravi- 
tational constant, and c is the velocity 
of light. To integrate the TOV equa- 
tions, one needs to know the equation 
of state for the entire expected density 
range of the neutron star, starting from 
the high density at the center to the 
surface densities. Therefore, we con- 
struct a composite equation of state for 
the entire neutron star density span by 
joining our equation of state for high 
density neutron rich matter to those 
with (i) 10^^ to 5 X 10^° g cm^^^ 16) ^j-^ 
5 X 10^° g cm~^ to lO"^ g cm"^,^'') and 
(iii) less than 10^ g cm~^. Thus we in- 
tegrate the TOV equations for the newly 
constructed equation of state and given 
central density e(r = 0) = €c with the 
boundary condition m(r = 0) = to 
give R and M. The radius R is defined 
by the point where P ~ 0, or, equiv- 
alently, e = eg, where eg (7.8g cm~^) 
is the density expected at the star sur- 
face. The total mass is then given by 
M = m{R). 

The results for the star structure parameters are listed in Table I and displayed 
in Fig. 3. This figure plots the mass as a function of the central density. The models 
are the same as in Fig. 2. 

In Table I, we also list additional parameters of interest. These are the moment 
of inertia 7, and the surface redshift z = , ^ = — 1 as a function of the central 

' ^/l-2GM/Rc^ 
density of the star. (For details, sec Rcf. 9)) These are important for the dynam- 
ics and transport properties of pulsars. From Fig. 3 and Table I, we see that the 
maximum masses of the stable neutron stars are 2.IM0, 2.2M0 and 2.5Mq and cor- 
responding radii are 12.1 km, 11.6 km and 13.8 km for MCH, NL3 and CH equation 
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3. The neutron star mass as a function 
radius. The models are the same as in 
Fig. 1. 
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Table I. Neutron star structure parameters. 
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of states, respectively. The corresponding central densities are 2.0 x 10^^ g cm~^ (> 7 
times nuclear matter density, where the nuclear matter density is 2.8 x 10^^ g cm^''), 
2.0 x 10^^ g cm^"^ (> 7 times nuclear matter density) and 1.5 x 10^^ g cm~^ (> 5 
times nuclear matter density) for MCH, NL3 and CH, respectively, at the maximum 
neutron star masses. Prom the neutron star structure point of view, our present 
model MCH is comparable to that of the NL3 model, which was constructed using 
recent heavy-ion collision data. The maximum masses calculated using our models 
are in the range of recent observations, " where the observational consequences 
are given below. 

Very recently, it has been observed that the best determined neutron star masses 
are found in binary pulsars, and they all lie in the range 1.35ibO.O4M0, except for 
those of the non-relativistic pulsars PSR J1012+5307, for which M = (2.1±0.8)Mq. 

There are several measured X-ray binary masses, and the heaviest among them 
arc Vela X-1 with M = (1.9±O.2)M0 ^o) and Cygnus X-2 with M = (1.8±O.4)M0. ^i) 
From the recent discovery of high-frequency brightness oscillations in low-mass X-ray 
binaries, the large masses of the neutron stars are around M = 2.3M0, M = 2.IM0 
and M = 1.9M0 in QPO 4U 1820-30, QPO 4U 1608-52 and QPO 4U 1636-536, 
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respectively. ^^-^ This provides a new method to determine the masses and radii of 
the neutron stars. Our results lie in the range of those predicted by the observational 
limits. -22) 

§4. Summary and conclusion 

We have discussed the SU{2) chiral sigma model, in which the isoscalar vector 
field is generated dynamically. In the present calculation, we have modified the 
SU{2) chiral sigma model by including two extra terms in the Lagrangian to ensure 
an appropriate value of the incompressibility at the saturation density of symmetric 
nuclear matter. By employing these and using a mean field approximation, we have 
obtained the nuclear equation of state at high densities. This is compatible with the 
equation of state that was derived using recent heavy-ion collision flow data. " 
The nuclear equation of state is reasonably softer than the original SU{2) chiral sigma 
model wc considered previously. We then employed the same nuclear equation of 
state in the neutron star matter calculation, where the composition of the star matter 
consists of neutrons, protons and electrons. The mass and radius of the neutron star 
were calculated from the TOV equations, and were found to be 2.1Mq and 12.1 km, 
respectively. The corresponding central density is 2.0 x 10^^ g cm~^ (> 7 times the 
nuclear matter density, 2.8 x 10^^ g cm~^) for the modified chiral sigma model. These 
values are comparable to the values calculated using the recent NL3 model, which 
was recently constructed using heavy-ion collision data. The overall structural values 
arc in good agreement with those given by the quasi-periodic oscillation method, 
a new method to determine the masses and radii of neutron stars. 

In the future we plan to investigate more systematically the properties of neutron 
stars in the presence of other baryons, i.e., hyperons and mesons with (hidden- 
)strangeness in the SU{3) chiral sigma model with the proper interactions between 
them derived from recent experiments. Moreover, we are interested in searching for 
a phase transition to quark matter (a new state of matter) by implementing this 
equation of state in heavy-ion collision simulation calculations. 
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